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ABSTRACT 

A  Galerkin  method  for  systems  of  PDE’s  in  circular  geometries  is  presented  with  motivating 
problems  being  drawn  from  structural,  acoustic  and  structural  acoustic  applications.  Depend¬ 
ing  upon  the  application  under  consideration,  piecewise  splines  or  Legendre  polynomials  are 
used  when  approximating  the  system  dynamics  with  modifications  included  to  incorporate 
the  analytic  solution  decay  near  the  coordinate  singularity.  This  provides  an  efficient  method 
which  retains  its  accuracy  throughout  the  circular  domain  without  degradation  at  the  singu¬ 
larity.  Because  the  problems  under  consideration  are  linear  or  weakly  nonlinear  with  constant 
or  piecewise  constant  coefficients,  transform  methods  for  the  problems  are  not  investigated. 
While  the  specific  method  is  developed  for  the  2-D  wave  equation  on  a  circular  domain  and  the 
equation  of  transverse  motion  for  a  thin  circular  plate,  examples  demonstrating  the  extension  of 
the  techniques  to  a  fully  coupled  structural  acoustic  system  are  used  to  illustrate  the  flexibility 
of  the  method  when  approximating  the  dynamics  of  more  complex  systems. 
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1  Introduction 


A  coniinoii  step  in  many  stnu  tnral.  aeonstie  ainl  fluid  applications  involves  the  solution  of  par¬ 
tial  differential  erpiations  (PDF/s)  modeling  the  physics  of  a  system  having  a  circidar  geometry. 
Exami)les  of  ap|)lications  in  which  this  is  important  iin  lnde  the  determination  and  control  of 
circular  plate  and  cylindrical  shell  dynami<s.  the  modeling  and  attenuation  of  noise  within  a 
cylindrical  cavity,  the  siimdation  of  Huid  flow  and  lumndarv  layer  growth  in  a  pijje,  as  well 
as  innumerable  other  applications  involving  circular  domains.  Due  to  the  complexity  of  the 
problems,  however,  analytic  solutions  usually  cannot  i)e  obtained  and  one  must  numerically 
approximate  the'  dynamics  under  consideration. 

From  separatieeii  of  variables  and  truncation  e)f  the  re'snlting  infinite  series,  aj)proximations 
in  the  form  of  modal  ('X])ansions  involving  Be'ssel  eomponents  can  occasionally  be  usenl  to 
approximate  PDF  dynamics.  However,  modal  expansions  can  be  determine'd  for  only  a  limited 
number  of  applications  aiul  hence  one  must  ofte;n  emi)loy  more  general  technique's  such  as  finite 
differene  e,  finite  element,  spline',  or  spe'ctral  expansions  in  order  to  approximate  the  dynatnie's 
of  the  system  uneler  consideration. 

In  the  cate'gory  of  spectral  me'the)els  for  circular  geometries,  (lalerkin.  e'ollocatieui  and  tan 
methods  have  be^en  stvufied  with  the  e  hoie'e  of  methoel  depending  upon  the  problem  being  e:ori- 
sidered.  Te)  date,  much  e)f  this  researe  h  has  eentered  arounel  the  simulation  e)f  fluiel  How  and 
boumiary  layer  growth  and  in  these  cases,  emphasis  has  usually  been  placed  on  collocation  due 
its  success  in  handling  complex  boundary  conditions,  variable  coefficients,  and  strotig  notdin- 
earities  [Ki.  17],  (lalerkin  methods  for  flows  on  spheres  are  discussed  in  [15],  but  this  is  done 
primarily  in  the  context  of  Fourier  expansions  involving  the  strong  form  of  the  modeling  flow 
equations  with  only  a  brief  discussion  concerning  Legendre  bases  being  included.  As  noted  in 
this  latter  reference,  care  must  be  exercised  when  a]>plying  any  of  the  previously  mentioned 
techniques  ((Jalerkin,  collocation  or  tan)  to  problems  with  coordinate  singularities  since  the 
incorrect  application  of  pole  conditions  can  significantly  degrade  the  accuracy  of  the  method 
as  well  as  introduce  strong  numerical  instabilities.  This  reference  also  includes  a  general  com¬ 
parison  between  the  residts  obtained  with  surface  harmonics  (eigenfunctions  of  the  Laplacian), 
modified  Robert  functions  and  Fourier  series  using  collocation  and  Calerkin  methods  in  the 
presence  of  the  coordinate  singularity.  The  use  of  the  mollified  Robert  functions  and  tech¬ 
niques  for  improving  their  conditioning  and  employing  fast  transforms  is  further  addressed  in 
[lOj.  Finally,  general  overviews  concerning  the  application  of  spectral  methods  to  problems 
with  coordinate  singularities  can  be  found  in  [1 1,  12].  We  point  out  that  the  emphasis  in  these 
references  centers  primarily  around  the  development  of  numerical  methods  for  flows  of  various 
types  and  hence  the  techniques  are  often  geared  toward  the  approximation  of  strongly  nonlinear 
operators  which  requires  transform  methods  to  facilitate  efficient  implementation.  While  these 
techniques  will  clearly  work  for  linear  problems  involving  circular  domains,  they  often  are  not 
optimal  nor  the  best  choice  for  the  applications  uniler  consideration. 

In  this  paper,  we  present  a  (Jalerkin  method  for  linear  or  weakly  nonlinear  problems  having 
circular  or  cylindrical  domains  using  piecewise  spline  and  spectral  bases.  Two  areas  from  which 
we  will  draw  examples  are  structural  dynamics  and  acoustics  (a  search  of  the  literature  reveals 
few  numerical  techniques  for  these  tyjies  of  problems  on  circular  domains).  Linearization  in 
the  first  case  is  often  justified  when  dealing  with  small  amplitude  vibrations  while  linearized 


acoustic  equations  are  often  employed  wlien  considering  acoustic  fields  having  sound  pressure 
levels  less  than  150  dB  (which  is  the  case  in  a  large  nuinher  of  acoustic  and  struct  ural  a(  (jiistu 
applications).  While  weak  nonlinearities  of  the  type  considered  in  [7]  can  also  he  ellicientls 
implemented,  the  discussion  here  will  concentrate  on  the  linear  ca.se.  Moreover.  th<‘  coeflici<'nts 
in  the  problems  of  interest  will  be  taken  to  be  constant  or  piecewise  constant  in  nature  (this  is 
certainly  a  reasonable  assumption  for  many  acoustic  problems  and  a  valid  comlition  in  a  large 
class  of  structural  applications).  Hence  we  will  not  address  the  use  of  transform  methods  for 
evaluating  variable  coefficients  and  noidinear  components. 

Throughout  the  discussion,  the  modeling  ecpiations  will  be  approximated  in  the  weak  or 
variational  form  in  accordance  with  energy  formulations  of  the  problems.  In  the  structural 
application,  this  is  done  to  reduce  smoothness  requirements  on  approximating  elements  and 
to  accommodate  structural  and  mateTial  discontinuities  as  well  as  unbounded  (discontinuous) 
input  operators.  The  consideration  of  the  acoustic  problem  in  weak  form  i)roves  to  lie  useful 
when  considering  coupled  structural  acoustic  .systems.  Moreover,  in  weak  form,  some  physit  al 
boundary  conditions  (e.g.,  the  hard  wall  conditions  which  we  consider  here)  are  untural  which 
implies  that  no  basis  alterations  are  necessary  when  implementing  the  method.  We  point  out. 
however,  that  if  the  application  warrants,  approximation  of  the  acoustic  dynattiics  using  the 
strong  form  of  the  system  equations  can  just  as  efficiently  be  accomplished  using  the  technicpies 
of  the  paper.  Moreover,  these  techniques  can  be  easily  modified  to  accommodate  a  tan  method 
if  dealing  with  more  complex  essential  boundary  conditions  (boundary  conditions  which  must 
be  explicitly  satisfied  in  a  variational  formulation  of  the  problem). 

In  the  second  section  of  this  paper,  a  2-D  acoustic  problem  on  a  circular  domain  is  con¬ 
sidered  with  a  Fourier- Legend  re  basks  being  u.sed  in  the  ensuing  approximations.  In  addition 
to  its  accuracy,  the  use  of  this  basis  proves  to  be  advantageous  in  control  applications  since  it 
facilitates  the  maintenance  of  uniform  stability  margins  under  approximation  [4j.  Moreov(‘r. 
the  the  use  of  the  spectral  basis  facilitates  approximation  in  a  quotient  space  which  is  often  the 
natural  state  space  in  acoustic  problems. 

The  techniques  are  then  applied  to  a  structural  application  in  Section  3  with  the  problem 
of  approximating  circular  plate  dynamics  being  used  to  motivate  the  analysis.  In  this  case, 
cubic  splines  modified  to  satisfy  boundary  conditions  are  used  in  conjunction  with  a  Fourier 
basis  to  obtain  a  suitable  finite-dimensional  approximation  to  the  problem.  The  use  of  splines 
rather  than  a  spectral  basis  is  motivated  by  smoothness  requirements  as  well  as  the  ease  with 
which  they  can  be  adapted  to  satisfy  essential  boundary  conditions,  and  can  easily  be  changed 
if  warranted  by  the  form  of  the  modeling  equations  and  boundary  conditions.  We  emphasize 
that  when  developing  a  general  technic|ue  that  could  be  applied  in  both  cases,  care  was  taken  to 
treat  the  coordinate  singularity  in  a  manner  that  prevented  degradation  of  accuracy  using  either 
basis  (spectral  or  piecewise  splines)  as  well  as  avoided  the  introduction  of  spurious  numerical 
instabilities. 

The  techniques  from  Sections  2  and  3  are  then  combined  in  Section  4  where  the  problem  of 
approximating  the  dynamics  of  a  fully  coupled  structural  acoustic  system  is  considered.  This 
demonstrates  the  flexibility  of  the  method  when  studying  a  more  complex  coupled  physical 
system  and  illustrates  the  manner  through  which  the  basic  techniques  can  be  extended  to  more 
general  problems.  The  convergence  properties  of  the  method  and  its  use  in  determining  physical 
properties  of  the  system  are  illustrated  through  a  set  of  examples.  Finally,  advantages  of  the 
method  in  control  applications  are  di.scussed  in  the  concluding  remarks. 
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2  The  Wave  Equation 


In  this  section  we  consider  the  ‘2-D  wave  equation  on  a  circular  domain  il  of  radius  a  with  the 
boundary  denoted  by  f.  For  Neumann  bouinlary  constraints  and  initial  conditions  (t)^  and  (pi, 
the  equations  of  motion  are 

<f)tt  =  c^A<f)  +  ,  (7-, 6^)  6  f?  >  0  , 

V(2!)-7‘  =  0  ,  (7-,  0)  6  r  ,  /  >  0  , 

(2.1) 

=  Mr.O)  ,  (r,0)Gn, 

=  <j,Ar,e)  ,  (i\e)en 

where  (j),  in  acoustics  applications  [14],  denotes  the  velocity  potential  and  c  is  the  speed  of 
sound.  We  note  that  the  Neumann  boundary  conditions  are  chosen  so  as  to  be  consistent  with 
the  hard  wall  conditions  in  the  physical  application  described  in  Section  4  and  can  be  altered 
to  fit  the  individual  problems  under  investigation. 

To  gain  insight  into  the  behavior  of  the  solution  and  to  motivate  the  form  of  the  approx¬ 
imation,  we  first  solve  the  equations  analytically  through  a  standard  separatioti  of  variables. 
For  the  homogeneous  problem  (g  =  0),  we  take  ^(f,7-, 0)  =  T(t)^(i\0)  to  arrive  at  the  two 
dimensional  Helmholz  equation 

7^4>  =  0  ,  ir,0)  €  , 

(2.2) 

V<l>-7-  =  0  ,(r,6>)Gr 

and  the  relation  T"  -f  =  0  ,  t  >  0  .  The  separation  constant  here  is  —7^  =  — (^)^  where  w 
is  the  circular  frequency  with  units  of  radians/sec.  To  find  4>,  we  separate  variables  once  more. 
Letting  4>(7‘,  0)  =  R{r)Q{0),  the  expansion  of  Helmholz’s  equation  yields  the  expression 

I  d  {  dR\  1  <{^e  ,  ^ 

vRdr  V'd7-  j  7-20 


which  implies  that  R  and  0  must  satisfy  the  differential  equations 

R'ia)  =  0  , 

and 

0"  -f  771^0  =  0 

0(0)  =  e(27r)  . 


R"+^-R'+  (7" 


(2.:i) 

(2.4) 


The  general  solution  to  (2.4)  is  periodic  and  is  given  by  Q[0)  =  m  =  0,  ±1,  ±2,  •  •  •  . 

Moreover,  the  general  solution  of  the  Bessel  equation  (2.;f)  is  R{r)  =  AiJm{v)  +  A^Aviiv') 
where  J,n  and  Y,n  are  the  nd''  Bessel  functions  of  the  first  and  second  kind,  respectively.  As 
noted  in  [1],  these  functions  display  the  asymptotic  bc'havior 


2”‘(77?.  -  1)! 


3 


for  fixed  m  and  '  — >  0.  Bec  ause  V’„  grows  unliouiidedlv  at  llu*  origin,  w<‘  lake  I  j  =  U  in  ordc'i 
to  guarantee  that  the  solution  rcunains  houndc'd  al  r  =  0.  I  he  dec  ay  ))ro|)cMl ic-s  of  will  l>e 
incorporated  when  developing  an  a|)proximation  mc-lliod  for  the  |)rol)lc‘m. 

Tlie  eigenvalues  of  (2.2)  are  them  clc-tc-rinined  hy  applying  the  Icoundary  eondilion  iunl 
solving  for  the  zeros  of  the  nonlin<-ar  e<|uatioti 


dr 


=  0  . 


(2.(i) 


The  corresponding  eigenvectors  or  mode's  are*  'I' h>f  "i  =  0-  ±1.  ±2.  •  •  ■. 
11  =  l,2,d,  •  •  •  which  implies  that  standing  wave's  have'  t  he'  form 


(Iji  .Sm ( 0  T  ^  ^  d f,^J^1  )( 


A  variation  of  parameters  apprejach  can  the'u  he'  use'cl  to  c'xtc'inl  this  to  t  he-  nonhomoge-ne'ous 
case.  With  suitable  smoothness  assumptieuis  on  the  initial  conditions  and  lore  ing  tc'rm.  it  c  an 
be  shown  that  the  .solution  to  (2.1)  has  the  form 


E  E  ^mti 


zrnO 


(2.7) 


m  =— 'X-  7i=:l 


where  the  c'oefficients  (i>,nn{f)  depend  on  nnel  the  time'-depende'ut  Fouric'r  coeffic  ic'uts  of 

g  (we  refer  the  reader  to  any  standard  PDE  te»xt  for  dc'tails  concerning  the'  juslilication  of  the 
formal  ralcuJation.s  outlined  here). 

Having  outlined  the  arguments  leading  to  the  analytic-  form  of  the  solution,  we'  now  want  to 
use  those  results  to  motivate  an  appropriate  CJalerkin  technie|ue  for  apjiroximating  the  .solution. 
To  facilitate  the  u.se  of  this  scheme  in  e-oupled  structural  acoustic  systems  (see',  for  e'xample, 
the  system  discussed  in  Section  4),  the  weak  form  of  the  .system  etpiations  will  he  c  onsidere'd. 
Moreover,  to  illustrate  the  use  of  the  metiuid  when  approximating  the  acoustic  potential  (whic  h 
is  determined  to  within  only  a  constant  [14]),  a  basis  suitable  for  a  cpiotient  space  will  be'  chosen. 

To  pose  the  system  (2.1)  in  weak  form,  we  take  the  state  to  be  0  in  the  siiace  H  =  L^(il) 
where  L^{fl)  is  the  quotient  space  of  over  the  constant  functions.  The  space  of  te'st  func  tions 
is  taken  to  be  the  quotient  space  V  =  The  inner  |)roducts  for  the  two  space's  are  take'u 

to  be 

(0,Ow  =  /  ^  V  =  f  V0-V^c/u; 

Ju  Ju 

(the  overbars  here  denote  comph'X  conjugation  as  compareil  with  those  usc'd  above*  to  de-lineate 

the  quotient  spaces).  Energy  cejusideratieuis  or  integration  by  parts  then  yields  the  variational 

form  _  _ 

I  +  /  V0-V^c/u;=  /  <i(,daJ  (2.S) 

Ju  Ju  j 

for  all  ^  in  l^. 

To  discretize,  we  begin  with  a  Fourie'r  expansieni  in  0  which  yie-lds  the*  approximate*  solution 


M 


Af 


iin$ 


(2.0) 
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We  point  out  that  the  use  of  the  complex  Fourier  expansion  simplifies  the  following  liiscussion 
both  in  describing  the  form  of  the  approximate  solution  and  the  construction  of  the  system 
matrices.  However,  when  combining  the.se  wave  results  with  those  of  the  circular  plate  to 
yield  an  approximation  scheme  for  the  coupled  system  (see  Section  4),  it  is  easier  to  use  a 
real  Fourier  expansion  when  performing  the  actual  computations  (this  is  due  to  the  presence 
of  the  piezoceramic  patches  on  the  plate).  The  interchange  between  the  two  expansions  is 
straightforward,  and  hence  details  concerning  the  implementation  of  the  real  Fourier  scheme 
are  left  to  the  reader. 

Several  possibilities  exist  for  spline,  Legendre  or  (Hiebyshev  expansions  of  </>„i(Lr)  both  in 
a  collocation  and  Galerkin  setting.  These  include  direct  expansions  which  maintain  the  parity 
of  the  solution,  refinements  to  incorporate  the  decay  of  the  .solution  at  the  origin,  and  mapped 
expansions  which  use  all  the  polynomials  and  yield  better  center  reso'ution.  These  expansions 
must  satisfy  the  condition 


d<f)^ 

~W 


=  0 


at  the  origin  which  guarantees  the  uniqueness  of  the  solution.  This  yields  the  requirement 


=  0  at  7-  =  0  ,  rn  ^  0  . 


(2.10) 


To  guarantee  differentiability  at  the  origin,  it  is  appropriate  to  require  that  the  remaining 
component  satisfies 


at  /•  =  0  . 


(2.11) 


As  detailed  in  [6],  one  expansion  of  which  satisfies  these  properties  is 


N” 


n=0 


where 


and 


m  — 


m 


jrn|  =  0,---,5 
|nj|  =  6,  •  •  • ,  M 


C(^) 


A(r)-  1/3 
Pnir) 


,  777  =  0,  77  =  1 

,  otherwise  . 


(2.12) 


Here  Pn{v)  is  the  77*^'  Legendre  polynomial  which  has  been  mapped  to  the  interval  (0,o).  The 
term  /’i(7‘)  —  1/3  when  in  =  0,77  =  1  results  from  the  orthogonality  properties  of  the  Legendre 
polynomials  and  arises  when  enforcing  the  condition  jQ<j)^[t^i\9,z)(Lj  =  0  so  as  to  guarantee 
that  the  functions  are  suitable  as  a  basis  for  the  quotient  space.  The  inclusion  of  the  weight 
7’l”‘l  incorporates  the  decay  of  the  analytic  solution  near  the  origin  (see  (2.5))  while  ensuring 
its  uniqueness  at  that  point.  Finally,  we  note  that  if  N  Legendre  functions  are  used,  the  limit 
is  given  by  N'"'  =  +  1  when  jrrr]  ^  0  and  =  N  when  m  =  0,  which  implies  that 

M  =  (2A/  +  1)(A^”*  +  1)  —  1  basis  functions  are  used  in  the  wave  expansion. 

Summarizing,  the  approximate  solution  to  (2.8)  can  be  expressed  as 


M  N”' 

</7^(f,r,0)=  Y.  E 

m=—M  n=I 


M 


£4(i)ef  (>■,») 


5 

i 

1 


(2.13) 


where  B^{r,6)  =  In  comparing  the  form  of  this  approximate  solution  with  the 

analytic  wave  solution  defined  in  (2.7),  it  can  be  seen  that  we  have  essentially  repla<  e<l  th<*  B<*s.sel 
components  by  weighted  Legendre  functions.  This  provides  an  approximate  .solutioti  which  i.> 
more  flexible  and  is  quicker  to  converge  in  a  variety  of  applications  than  that  obtained  with  a 
purely  eigenfunction  expansion  (indeed,  for  more  complex  systems  involving  wave  components, 
the  eigenfunctions  are  unavailable  and  would  have  to  be  apjjroximated  before  they  could  be 
used  as  a  basis). 

To  provide  an  approximating  finite-dimensional  space  for  the  problem,  we  define  the  sub¬ 
space  =  span{5if^}  and  take  the  product  space  for  the  first-order  probleiTi  to  be  = 
X  The  restriction  of  the  infinite-dimensional  system  (2.8)  to  the  space  x 

then  yields 

/  I  [  -cj^du> 

Jii  Jq  Jq 

for  ^  in  .  The  corresponding  matrix  system  is 

d^{t)  ]  r  0  ' 

0  d^{t)  ^  G^{t) 

(2.14) 


where  t?'^(Z)  =  ^2(^)1  •  •  •  ?  denotes  the  A4  x  1  vector  containing  the  approximate 

state  coefficients.  The  component  matrices  and  vectors  are  given  by 

[Gf  (Oj,  = /^  ,  (2.15) 

=  /„  ■  W'iu.  =  (^,  , 

where  the  index  ranges  are  k,£  =  We  note  that  the  presence  of  the  matrix  in 

the  1?"'^  component  of  (2.14)  and  mass  matrix  is  due  to  the  form  of  the  (/-inner  product. 

We  point  out  that  the  Fourier-Galerkin  technique  leading  to  the  matrix  system  (2.14)  is  quite 
general  in  nature  and  can  be  used  to  approximate  wave  solutions  in  cases  where  the  analytic 
solutions  described  at  the  beginning  of  this  section  are  unavailable  or  impractical  to  use.  The 
general  techniques  described  here  are  also  easily  adapted  to  different  boundary  conditions  and 
more  complex  systems  involving  acoustic  components.  Test  cases  illustrating  the  basic  method 
are  given  in  the  following  examples  and  the  extension  of  these  techniques  to  a  more  complex 
coupled  structural  acoustic  system  is  presented  in  Section  4. 


■  0  d^(t)  r  0 

0  \[  d^{t)  J  ~  [  -K^ 

'  0  1  r  i?^(o)  1  _  [  5)^  ■ 

0  t9^(0)  ~ 
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2.1  Example  1 

As  a  first  test  of  this  approximation  metliod,  we  considered  the  forced  wave  equation 

</)((  =  c^A<c>  + 5(<,r,0)  ,  (r,9)eilj>0. 

V(t>-f  =  0  ,  (r,0)6r.<>O. 

4>{0,rJ)  =  (l>t(0,rj)  =  0  ,  (r,0)€n 

on  the  circular  domain  of  radius  «  =  1.  The  true  solution  r,  0)  =  <^[cos(27r7  )—  l]^[sin(d0)  + 
cos(0)]  was  used  to  generate  the  forcing  function  g,  and  the  problem  was  then  discretized  to  yield 
the  matrix  system  (2.14).  The  absolute  and  relative  errors  obtained  by  integrating  this  .system 
to  time  T  =  1  and  comparing  with  the  known  true  solution  on  a  uniform  mesh  are  recorded 
in  Table  I.  The  true  solution  at  that  time  is  depicted  in  Figure  1  while  line  plots  comparing 
the  approximate  and  true  solutions  along  the  central  line  =  {(x,?/):  —1  <x<  l,t/  =  0}  are 
given  in  Figure  2. 

The  solution  for  this  problem  was  chosen  in  a  manner  that  allowed  us  to  fix  the  Fourier 
limit  and  examine  the  convergence  as  the  CJalerkin  limit  was  increased.  The  rapid  convergence 
exhibited  by  the  results  in  Table  1  is  consistent  with  that  expected  from  the  radial  Legendre 
basis.  A  slight  decrease  in  the  convergence  rate  is  noted  when  N  =  15,  and  we  believe  that  this 
is  due  to  error  in  time  discretization  rather  than  spatial  approximation  (a  standard  fourth/fifth 
order  Runge-Kutta  routine  was  used  to  integrate  the  system).  In  this  and  other  examples 
that  we  have  examined,  it  can  be  noted  that  as  long  as  the  solution  being  approximated  is 
sufficiently  smooth  at  the  origin,  the  approximation  method  performs  well  and  accurate  results 
can  be  obtained  with  a  relatively  small  number  of  basis  functions.  For  example,  the  plots  in 
Figure  2  demonstrate  that  while  small  oscillations  near  the  origin  are  present  in  the  approximate 
solutions  obtained  with  N  —  ‘i  and  N  =  6,  they  are  gone  with  larger  discretization  limits  and 
if  fact,  the  approximate  solution  obtained  with  Af  =  9  is  graphically  indistinguishable  from 
the  true  solution.  For  less  smooth  functions  (e.g.,  functions  that  are  only  continuous  and 
differentiable  at  the  origin),  some  loss  of  accuracy  does  occur,  although,  in  applications  this 
has  not  been  a  factor  since  these  types  of  discontinuities  generally  do  not  occur  in  the  physical 
systems  whose  dynamics  we  are  simulating. 


M 

N 

2M 

ll^irue  ^appll 

3 

3 

54 

1.0067-0 

.1357-0 

3 

6 

96 

.1453-0 

.1959-  1 

3 

9 

138 

.1269-2 

.1710-3 

3 

12 

180 

.8.567  -  4 

.11.55-4 

3 

15 

222 

.1912-4 

.2578  -  5 

Table  1.  Absolute  and  relative  errors  when  approximating  the  wave  solution 
(f>{t,r,d)  =  f^[cos(27rr)  —  l]^[sin(30)  +  cos{0)]. 
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2.2  Example  2 

A  second  means  of  testing  the  accuracy  of  the  iiitiss  and  stiffness  matrices  in  the  system  (2.14) 
and  hence  checking  the  accuracy  and  convergence  properties  of  the  approximation  method  is 
by  discretizing  the  eigenvalue  problem 

+  7^^  =  0  ,  G  n  . 

V<Hr==0  ,(r,<?)Gr 

(see  (2.2))  which  arises  when  separating  variables  in  the  homogeneous  wave  ecjuation.  .-Xs 
discussed  previously,  the  eigenvalues  7„4„  are  determined  by  solving  for  the  zeros  of  the  nonlinear 
system  —  0  where  J„i  is  the  Bessel  function  of  the  first  kind.  Several  values  obtained 

with  a  —  I  are  summarized  in  Table  2  (see  page  .'M.'f  of  [9]  for  details). 

To  compare  these  results  with  those  obtained  via  the  Fourier-CIalerkin  expansion,  it  is  noted 
that  under  approximation  in  weak  form,  the  eigenvalue  problem  (2.2)  yields  the  generalized 
matrix  eigenvalue  problem 

(2.16) 

with  =  «  =  1  in  the  mass  and  stiffness  matrices  and  respectively  (see  (2.15)  for  the 
definition  of  these  matrices).  The  approximate  eigenvalues  obtained  by  solving  (2.16)  are  sum¬ 
marized  in  Tables  3  and  4  for  the  limit  choices  M  =  6,  iV  =  6  and  M  =  6,  A  =  12,  respectively. 
With  the  first  choice  of  limits,  it  becomes  difficult  to  distinguish  the  higher  order  eigenvalues 
and  these  are  omitted  from  the  table.  We  first  note  that  with  M  =  6,  =  12,  very  accurate  ap¬ 

proximations  are  obtained  with  the  largest  relative  error  (4.4  x  10"'’)  occurring  when  m  =  n  =  6. 
Moreover,  we  see  that  while  the  definition  of  the  modified  Fourier  component  m  (see  (2.12)), 
which  incorporates  the  analytic  decay  of  the  Laplacian  near  the  origin,  changes  at  rn  =  5,  this 
has  not  reduced  the  accuracy  of  the  method.  Hence  this  example  further  illustrates  the  efficiency 
and  accuracy  of  this  approximation  method  for  the  wave  equation  on  a  circular  domain. 


n 

VI  =  0 

VI  =  1 

VI  =  2 

m  —  3 

VI  =  4 

VI  =  5 

VI  =  6 

0 

1.8412 

3.0542 

4.2012 

5.3176 

6.4156 

7.5013 

1 

3.8317 

5.3314 

6.7061 

8.0152 

9.2824 

10.5199 

11.7349 

2 

7.0156 

8.5363 

9.9695 

11.3459 

12.6819 

13.9872 

15.2682 

3 

10.1735 

11.7060 

13.1704 

14.5858 

15.9641 

17.3128 

18.6374 

4 

13.3237 

14.8636 

16.3475 

17.7887 

19.1960 

.20.5755 

21.9317 

5 

16.4706 

18.0155 

19.5129 

20.9725 

22.4010 

23.8036 

25.1839 

6 

19.6159 

21.1644 

22.6716 

24.1449 

25..5898 

27.0103 

28.4098 

Table  2.  Values  of  7„„i 

obtained  from  the  Bessel  condition 

dJmh)  _ 

dr 

n 

VI  =  0 

VI  =  1 

VI  —  2 

VI  =  3 

VI  =  4 

m  =  5 

VI  =  6 

0 

1.8412 

3.0542 

4.2012 

5.3176 

6.4156 

7.5013 

1 

3.8317 

5.3314 

6.7061 

8.0152 

9.2824 

10.5199 

11.7349 

2 

7.0159 

8.5378 

9.9714 

11.3472 

12.6823 

13.9871 

15.2709 

3 

10.2668 

11.7350 

13.1858 

14.6280 

16.0386 

17.4007 

18.7540 

4 

13.7669 

15.2709 

17.2673 

18.3866 

Table  3.  Values  of  7,„„  obtained  with  M  =  6,  )V  =  6  Fourier-Galerkin  basis  functions. 
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n 

m  =  0 

BBDI 

isea 

ni  =  3 

771  —  4 

Itl  =  0 

711  =  6 

0 

1.8412 

3.0542 

4.2012 

5.3176 

6.41.56 

7..5013 

1 

3.8317 

5.3314 

6.7061 

8.0152 

9.2824 

10.5199 

11.7.319 

2 

7.0156 

8.5363 

9.9695 

1 1 .3459 

12.6819 

13.9872 

15.26S2 

3 

10.1735 

11.7060 

13.1704 

14.58.58 

15.9641 

17.3128 

18.6374 

4 

13.3237 

14.8636 

16.3475 

17.7888 

19.1960 

20.57.55 

21.9317 

5 

16.4717 

18.0159 

19.5131 

20.9730 

22.-i018 

23.8043 

25.1860 

6 

19.6250 

21.1860 

22.6929 

24.1.5.59 

25.5937 

27.01.52 

2S.4223 

Table  4.  Values  of  7„„i  obtained  with  M  =  6,  =  12  Fourier-CIalerkiii  basis  fiiiictiuiis. 


3  The  Plate  Equation 

A  second  area  in  which  one  must  commonly  approximate  the  dynamics  of  PDE-based  models 
posed  on  circula”  domains  occurs  when  numerically  simulating  structural  dynamics.  While  the 
assumption  of  small  amplitude  vibrations  often  leads  to  linear  models,  these  models  are  often 
sufficiently  complex  so  as  to  warrant  special  numerical  techniques  to  capture  the  physics  of  the 
system.  As  a  motivating  example,  we  consider  a  thin  circular  plate  having  piezoceramic  patches 
bonded  in  pairs  to  its  surface  (see  Figure  3).  When  a  voltage  is  applied  to  the  patches,  stresses 
are  generated  which  can  be  used  to  invoke  in-plane  forces  and/or  bending  moments  in  the 
underlying  structure  [8].  In  this  manner,  the  patches  can  be  used  to  control  plate  vibrations  [2] 
or  acoustic  sound  pressure  levels  when  the  plate  is  an  active  component  in  a  structural  acoustic 
system  [3,  4,  5].  For  the  analysis  which  follows,  we  assume  that  the  edge  of  the  plate  is  clamped 
since  this  quite  closely  approximates  the  condition  found  in  several  applications  of  ititerest. 
We  point  out,  however,  that  the  techniques  which  follow  are  easily  adapted  to  other  essential 
boundary  conditions  if  warranted  by  the  physical  model. 


Piezoceramic  Patches 


Figure  3.  A  thin  circular  plate  with  piezoeceramic  patches  bonded  in  pairs  to  its  surface. 
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For  a  thin  circular  plate  of  radius  a  having  clamped  edges,  the  strong  form  of  the  equations 
modeling  the  transverse  motion  is 


if  Mr  2dMr  \  dMe  'lifMro  2  dMrs  1 


+ - : 


+ -g;5- +  “ST 

w{t,a,0]  =  ^(<,0,0)  =  0 
or 

w{Q,i\6)  =  Wq{i\0)  ,  wt{Q,i\d)  =  wx{r,d) 


+ 


r  di-dd  do 


+ 


=  f{t,r,0)  , 


(.11) 


where  pp  and  h  denote  the  density  and  thickness  of  the  plate  and  w  is  the  transverse  displace¬ 
ment.  The  general  moments  are  given  by 

.Ad  7*  —  A/f  Mr}pe 
Me  =  —  (Me)pe 

MrS  =  Mr0 


where 


d^w  u  dw  u  d^w  \ 


^  dr^  ^  r  dr  r'^  dO^  /  ^  ^\br‘^dt  ^  r  drdt  r'^  dO^dt ) 


dfw  ty  d^w  V  d^w 


I  dw  1  d^w 
Me  =  D\ -  — + 


r  dr  r^  d0^ 


d^w\  ( 1  d^w 


1  cfw  (Pw  \ 

r  drdt  '  7-2  d0'^dt  dr^dt  j 


(3.2) 


r  drd0 


1  dw\  ./l  d^w  1  d^w\ 

drd0dt 


are  the  internal  plate  moments  (including  discontinuous  changes  in  D,i/  and  co  due  to  the 
bonding  of  the  patches  to  the  plate  -  see  [8]),  and 

{Mr)pe  =  {Me)pe  =  ^  K:iUi{t)xi{r,  0) 

i=l 

are  the  applied  moments  generated  by  s  pairs  of  patches.  With  E  denoting  the  Young’s  modulus, 
the  parameters  j/,  D  =  represent  the  Poisson  ratio,  flexural  rigidity,  and  damping 

coefficient  for  the  plate/patch  structure.  Here  denotes  the  characteristic  function  which 

has  a  value  of  1  in  the  region  covered  by  the  patch  and  is  0  elsewhere.  Moreover,  Ui{t)  is  the 
voltage  into  the  patch  and  K'i  is  a  parameter  which  depends  on  the  geometry,  piezoceramic 
material  properties  and  piezoelectric  strain  constants  (see  [8]  for  details).  We  point  out  that  the 
piezoceramic  material  parameters  K'i,i  =  1,  •  •  ,s  as  well  as  the  plate  parameters  pp,  D,C[)  and 
should  be  considered  as  unknown  and  in  applications  must  be  estimated  using  data  fitting 
techniques  analogous  to  those  discussed  in  [7|.  The  piecewise  nature  of  the  material  parameters 
and  input  moments  is  one  motivation  for  approximating  the  problem  in  a  weak  or  variational 
form. 

To  motivate  the  development  of  a  Fourier-Galerkin  scheme  for  this  problem,  we  proceed  in 
a  manner  analogous  to  that  used  in  the  last  section  and  consider  first  a  simplified  version  of 
the  problem  where  separation  of  variables  can  be  employed.  To  this  end,  we  analyze  first  the 
equation 

Pph—  +  DV^w  =  0  .  (3.3) 
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which  models  the  transverse  motion  of  an  undamped  thin  circular  plate  devoid  of  patches  and 
having  constant  stiffness  and  density.  Because  the  manner  of  solution  is  similar  to  that  discussed 
in  the  last  section,  we  highlight  only  those  details  which  are  necessary  for  motivating  the  form 
of  the  approximate  solution  for  this  problem. 

As  before,  we  separate  the  temporal  and  spatial  variables,  in  this  ca.se  taking  w{t,r,0)  = 
T{t)W{r,d),  to  obtain 

dw  (•^•4) 

Wiaj)^—{a,e)  =  0 

and  T"  +  (jj^T  =  0  ,t  >  0.  The  separation  constant  is  7^  =  where  ui  again  denotes  the 

circular  frequency.  A  second  separation  of  variables  then  yields  solutions  of  the  form 


where  and  Ku  are  Bessel  functions  of  the  first  and  second  kind  and  /„i(-)  =  and 

-p{z)  —  Ip{z)]/  sm(p-K)  are  modified  Bessel  functions  of  the  first  and  second 
kinds,  respectively.  We  note  that  for  fixed  m  and  z  0,  these  functions  have  the  limiting 
forms 


~  2-7n! 


2”‘m! 

2”^(m-  I)! 
2z’'‘ 


(3.6) 


(see  [1])  which  immediately  shows  that  the  Bessel  function  and  modified  Bessel  function  of  the 
second  kind  grow  unboundedly  at  the  origin.  Hence  C,n  and  Dm  are  taken  to  be  zero  in  order  to 
avoid  infinite  deflections  and  stresses  at  r  =  0.  Furthermore,  the  condition  W{a,0)  =  0  yields 
the  condition 

/m(A) 


A  —  —A? 

—  ^7u  j  /  »  V 

1  ^  j 


(3.8) 


where  A  =  7a.  Finally,  the  boundary  condition  ^(a,0)  =  0  requires  that  A  solve 

/..(A)^^  -  ■'..(A)^^  =  0 

in  order  that  (3.5)  yields  a  nontrivial  solution  to  (3.4).  By  solving  (3.8),  one  obtains  the  eigen¬ 
values  A„i,i  and  by  combining  (3.5)  and  (3.7)  with  the  Fourier  coefficients,  the  eigenfunctions 
or  modes  are  found  to  be 


J„ 


Amur 


J m  (  Anin  ) 


(3.9) 


^  /  Au  (•^mn)  \  ^ 

With  suitable  smoothness  assumptions  on  the  the  initial  conditions,  it  then  can  be  shown  that 


the  solution  to  (3.3)  has  the  form 


00  00 


w 


{t^r,0)=  ^  5Iu;„m(f) 

m=:— 00  n=t 


J  ( 'A71  (^vnu)  j  ( 

'  1“^ }  ~  ImiKnn) 


(3.10) 


where  the  coefficients  Wmn(i)  depend  on  the  initial  conditions.  Again,  we  note  that  the  motiva¬ 
tion  given  here  is  formal  and  we  refer  the  reader  to  any  standard  PDF  tex  for  details  concerning 
this  derivation. 


Having  motivated  the  general  form  of  the  analytic  solution  for  a  specific  case  of  the  plate 
equation,  we  can  now  present  a  Fourier-dalerkin  method  for  approximating  plate  dynamics  in 
the  more  general  cases.  As  noted  at  the  beginning  of  the  section,  the  weak  form  of  the  system 
equation  (3.1)  will  be  used  to  accommodate  potentially  discontinuous  material  parameters. 
Moreover,  the  use  of  the  weak  form  accommodates  unbounded  input  operators  (as  is  the  case 
when  piezoceramic  patches  are  used  as  control  elements)  as  well  as  red\ices  the  smoothness 
requirements  on  the  basis  functions. 

The  state  for  the  problem  is  taken  to  be  w  and  the  Hilbert  space  H  —  f'^(ro)  with  the 
energy  inner  product 

{iiuv)h  =  /  Ppt^wTid-i 

y  Pq 

is  used  as  the  state  space.  We  also  define  the  Hilbert  space  of  test  functions  V  =  Hq[Tq)  = 
{%)  G  //^(Fo)  :  i){a,0)  =  V>^(rt,0)  =  0}  with  the  inner  product 


where  cp  =  0  in  the  moments  appearing  in  this  definition  and  (F,  =  fro  FCid'y  with 

d-y  =  rdrdO.  As  detailed  in  [6],  a  weak  or  variational  form  of  the  equation  describing  the  motion 
of  a  damped  thin  circular  plate  having  s  piezoceramic  patch  pairs  is 

PphwuVd^  +  + 

"'T 0  JTo 

for  all  T]  ^  V.  Again,  the  internal  moments  are  given  in  (3.2). 

To  obtain  an  appropriate  form  for  the  approximate  solution,  we  begin  with  the  Fourier 
expansion  in  6 

M 

w^{Fr,0)= 

m=—M 

As  was  the  case  when  describing  the  wave  approximation,  we  will  use  the  complex  Fourier 
expansion  while  describing  the  method  since  it  simplifies  the  notation.  Due  to  the  potential 
presence  of  patches  or  other  actuators  on  the  plate,  however,  the  method  is  more  easily  imple¬ 
mented  using  the  real  (trigonometric)  expansion  (with  patches  present,  the  complex  expansion 
leads  to  a  complex  system  matrix  which  proves  troublesome  when  solving  the  Riccati  equation 
in  the  control  problem).  We  have  omitted  details  concerning  the  real  expansions  used  in  the 
implementation  since  it  is  straightforward  to  interchange  between  the  two  expansions. 

To  determine  an  appropriate  expression  for  r),  we  first  note  that  it  must  satisfy  the 

conditions  ^  _ 

=  0  ,  m  7^  0  ;  =  0 

at  ?■  =  0  in  order  to  guarantee  uniqueness  and  differentiability  (see  the  discussion  about  the 
analogous  requirement  for  the  wave  solution).  In  light  of  these  requirements,  an  appropriate 
expansion  for  Wm(t,r)  is 

JSfm 

«’m(<,7’)  =  «’mH(0r’"‘’F,"‘(r) 

71=1 


L 


—M 


t9 


dff^  ^  drdo 


de 


dy 


(3.11) 
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where  B"‘(r)  is  the  modified  rtihir  spline  satisfying  =0  with  the  comliliun 

— =  0  being  enforced  when  ?n  =  0  (this  latter  <'unditio?i  gnarante(\s  dith^rentiabililv  at 
the  origin  and  implies  that  N''^  =  J\f  l  when  ?n  ^  0  and  A'”*  =  N  when  m  =  0,  where  N 
denotes  the  number  of  modified  eubir  splines).  The  total  number  of  jrlate  basis  funetioiis  is 
M  =  (2M  +  1)(A^  +  1)  —  1.  The  inclusion  of  the  weighting  term  with 


f  0  ,  =  0 

m  —  < 

I  1  ,  77*  /  0 

is  motivated  by  the  asymptotic  behavior  of  the  Bessel  functions  (see  (d.fi))  as  r  0.  It  also 
serves  to  ensure  the  uniqueness  of  the  solution  at  the  origin.  The  Fourier  coeHicieiit  in  the  weight 
is  truncated  to  control  the  conditioning  of  the  mass  and  stiffness  matrices  (see  e.xamples  in  [(>]). 
To  simplify  notation,  the  approximate  solution  is  written  as 


M  N'^ 

m=—M  n=l 


*:=] 


(.•1.12) 


where  Bjf{7\0)  =  7-l”*li9”‘(7')e‘”‘^.  A  comparison  with  (3.10)  indicates  that,  as  was  the  ca.se  with 
the  wave  approximation,  the  Bessel  components  of  the  analytic  solution  are  essentially  replaced 
by  approximating  elements  suitable  for  the  problem  under  consideration.  In  the  wave  problejn. 
the  approximate  solution  was  expressed  in  terms  of  weighted  Legendre  polynomials  whereas  here 
the  radial  basis  is  comprised  of  weighted  cubic  splines.  In  both  cases,  the  weights  were  used  to 
incorporate  the  analytic  solution  behavior  at  the  origin  into  the  approximate  solution,  thus  guar¬ 
anteeing  that  it  had  the  correct  continuity  and  differentiability  characteristics  at  that  point. 

To  obtain  a  corresponding  matrix  system,  the  M  dimensional  approximating  subspace  is 
taken  to  be  =  span{B^}  and  the  product  space  for  the  first-order  system  is  x  . 
The  restriction  of  the  infinite-dimensional  system  (3.11)  to  the  space  x  then  yields  the 
matrix  equation 


A'^ 

0 

0 

, 

1  1 

1  1 

1 

1 

r  0 

- 

j.'JV 

-K^ 

J 

■  w^it)  ■ 

vv^(0 

-f 

0 

Bf 

u{t)  -1- 

1  1 

1 _ » 

0  ■ 

■  >V^(0)  ■ 

’  9\  ' 

0 

VV^(O) 

.  92  . 

where  ■  ■  ■  ,  u;;v'(0]^  denotes  the  A/'x  1  vector  containing  the  approximate 

state  coefficients.  The  component  matrices  and  vectors  are  given  by 

A  ^  =  A  £ji  -|-  A  p2  +  A  D3  -|-  A  D4  -|-  A  05 

A^^j  A(,pj  -f-  I\fp2  T  Af.^3  Af./,4  T  A(.p5 
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where 


=  L  +  7-^  +  T^^J  ’ 

r,..  1  f  r.\i  dBf  ,  1  ff^Bf  ,  I'd^Bf]  dB^  , 
t^D2]#.fc  Jr,^[r^  dr  30^  ^  r  dr^  dr  ' 

[A'nl  =  /  nU^  +  l^  +  ii^l^d^ 

‘  Ju  ^  [r^  dr  ■*■  r-*  30^  r2  dr^  30^  ’ 

r..  1  _  o  /  r..,  X  f  I  1  dBf]  d^B?  , 

[^D4]t,k  ‘'^[ri  drd0  r^  30  \  drd0  ' 

fL^  1  o  /  r^/.  ^\  i  d^Bf  I  dBf]  dW  ^ 

[Ads],.,  =  2  D(  I  -  .)  |-^—  +  -_J  rf, 

with  the  index  ranges  are  k,i  =  1,  •  •  •  ,M.  The  matrices  Kcj,\  —  are  defined  similarly  with 
the  inclusion  of  the  parameter  co  in  the  various  integrals.  Finally,  we  remind  the  reader  that 
Pp,  D,  V  and  cq  are  piecewise  constant  in  these  definitions  due  to  the  presence  of  the  patches. 

For  application  purpose,  it  useful  to  note  that  the  matrix  system  for  the  plate  can  be  written 
as  the  Cauchy  system 

y^{t)  =  A^y^it)  +  B^u{t)  +  G^{t) 

where y^(<)  =  [W^(<), VV‘^(f)]  =  denotes  the ‘ZA/* x  1  vector 

containing  the  generalized  Fourier  coefficients  for  the  approximate  displacement  and  velocity. 
In  this  form,  finite-dimensional  parameter  estimation  and  control  problems  can  be  readily  dis¬ 
cussed. 


3.1  Example  3 

As  an  initial  test  of  the  Fourier- Galerk in  approximation  scheme,  we  considered  the  undamped 
steady  state  problem 

d^Mr  2dMr  IdMe  ‘Zd^M^e  2  dMre  1  d^Mp  _ 

dr^  r  dr  r  dr  **"  r  dr^  r^  dO  dO"^  ’ 

with  the  moments  given  by  (3.2).  In  these  expressions,  the  Poisson  ratio,  flexural  rigidity  and 
damping  coefficient  were  taken  to  be  =  |,  £)  =  |  and  cp  —  0,  respectively.  The  plate  radius 
was  taken  to  be  a  =  .6  and  the  true  solution  w{r,9)  =  (cos(27rr/a)  —  l)sin(^)  was  used  to 
generate  the  forcing  function 

f{r,6)  =  ^  ^3  -j-  (^—6irr/a  -|-  Wir^r^fa^'j  sm{2irr/a) 

-I-  (—3  -I-  IZir^i'^fa^  +  cos(27rr/a)l  sin(0)  . 
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The  absolute  and  relative  errors  for  this  example  are  summarized  in  Table  5.  By  choosing  a  true 
solution  with  a  known  Fourier  coefficient,  the  Fourier  limit  could  be  fixed  at  A/  =  1  and  the  con¬ 
vergence  rate  could  be  observed  as  the  Galerkin  limit  N  wcis  increased.  For  sufficiently  smooth 
forcing  functions,  we  would  expect  the  method  to  exhibit  0{h*)  convergence  as  a  result  of  the 
cubic  spline  basis  [18,  19].  To  check  this,  asymptotic  errors  were  calculated  by  dividing  the  pre¬ 
vious  relative  errors  by  16.  Since  the  number  of  radial  basis  functions  is  doubled  each  time,  this 
provides  a  means  of  checking  whether  or  not  the  expected  convergence  rate  is  being  maintained. 
By  comparing  the  results  in  the  and  columns  of  the  table,  see  that  0{h*)  convergence 
is  being  exhibited  by  the  method  thus  providing  an  initial  test  in  the  efficacy  of  the  method. 


M 

N 

size(A-^) 

Untrue  ^appll 

mi 

Asyni.  Error 

1 

5 

15  X  15 

.4362  -  2 

.2181-2 

1 

10 

30x  30 

.1454-2 

.7269  -  3 

.1363-3 

1 

20 

60x  60 

.6146  -  4 

.3073  -  4 

.4543  -  4 

1 

40 

120  X  120 

.6717-6 

.3358  -  6 

.1921  -  5 

Table  5.  Absolute  and  relative  errors  for  approximate  solutions  for  Example  1. 


3.2  Example  4 

A  second  means  of  testing  the  accuracy  and  efficiency  of  the  approximation  method  is  by  dis¬ 
cretizing  the  eigenvalue  problem  (3.4)  and  comparing  the  approximate  eigenvalues  and  eigen¬ 
functions  with  analytic  values  that  have  been  calculated  for  a  simplified  structure.  This  also 
provides  a  means  of  determining  analytic  values  of  the  natural  frequencies  of  a  uniform  plate 
which  can  then  be  compared  with  values  obtained  numerically  and  experimentally  for  structures 
involving  plates  to  which  patches  have  been  bonded  [2]  and  plates  that  have  been  incorporated 
in  structural  acoustic  systems  (as  described  in  Section  4). 

For  a  uniform  undamped  plate  to  which  no  patches  are  bonded,  approximation  of  the 
eigenvalue  problem  (3.4)  via  the  Fourier-Galerkin  method  described  in  this  section  yields  the 
generalized  matrix  eigenvalue  problem 

=  (3.15) 

with  the  mass  and  stiffness  matrices  and  Kp  defined  in  (3.13).  The  plate  dimensions 
a  =  .2286m  (9in),  pp  =  2700kgl'm^,  h  =  .00127m  (.05m),  and  parameter  choices  E  = 
7.1  X 10’'’  Nfm^  and  i/  =  .33  were  used  which  then  yields  the  flexural  rigidity  D  =  13.6007  N  •  m. 
We  point  out  that  these  choices  are  consistent  with  the  dimensions  and  parameters  of  an 
experimental  plate  currently  being  used  in  the  Acoustics  Division,  NASA  Langley  Research 
Center  so  that  the  frequencies  determined  here  could  be  compared  with  those  obtained  obtained 
experimentally  (again,  see  the  application  in  the  following  section). 

By  noting  the  relationship  /  =  where  /  is  the  frequency  expressed  in  hertz  and  uj  is 
the  circular  frequency,  the  natural  frequencies  of  the  fixed  circular  plate  can  be  written  as 
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For  the  given  dimensions  and  parameter  values,  several  frequencies  deriving  from  the  Bessel 
solutions  A  to  the  nonlinear  equation  (3.8),  as  reported  on  page  8  of  [13],  are  given  in  'Fable  (i. 
Approximate  frequencies  obtained  by  solving  the  generalized  matrix  eigenvalue  problem  (3.15) 
with  the  basis  limits  M  =  6  and  N  =  24  are  tabulated  in  Table  7.  We  point  out  that  in  these 
tables,  the  Fourier  number  m  can  be  interpreted  as  the  number  of  nodal  diameters  while  n  is 
the  number  of  nodal  circles,  not  including  the  boundary.  In  comparing  the  Bessel  and  Clalerkin 
results  in  the  two  tables,  it  can  be  seen  that  the  Fourier-Galerkin  provides  accurate  mass  and 
stiffness  matrices  which  translates  to  accurate  approximations  to  the  natural  frequencies  for 
the  plate.  With  the  basis  limits  used  here,  the  largest  relative  error  for  the  frequencies  shown 
here  is  8.1  x  10“'*  when  m  =  3,  n  =  6.  Moreover,  we  see  that  while  the  definition  of  m  in  the 
weighting  term  changes  at  m  =  1,  the  method  loses  no  accuracy  at  that  point. 

To  qualitatively  test  the  ability  of  the  Fourier-Galerkin  method  to  accurately  approximate 
the  decay  of  the  solution  in  neighborhoods  of  the  origin  for  increasing  m,  we  compared  the 
analytic  eigenfunctions  given  by  (3.9)  with  those  obtained  via  the  Fourier-Galerkin  method 
with  M  =  5  and  N  =  16.  Representative  results  using  the  two  techniques  are  plotted  in 
Figures  4  and  5  with  additional  examples  demonstrating  the  3-D  behavior,  corresponding  2-D 
slices  and  error  results  given  in  [6].  Qualitatively,  the  shape  of  the  eigenfunctions  in  the  two 
sets  can  be  seen  to  be  graphically  identical  and,  by  comparing  the  results  for  the  (rn  =  5,  n  =  0) 
mode,  it  can  be  seen  that  the  approximation  scheme  is  accurately  capturing  the  behavior  of  the 
solution  near  the  origin.  In  combination  with  the  eigenvalues  results  listed  in  the  tables,  this 
demonstrates  that  the  approximation  method  accurately  captures  the  physics  of  the  problem 
throughout  the  circular  domain  with  no  loss  of  accuracy  resulting  at  the  origin,  in  spite  of 
the  coordinate  singularity.  Moreover,  as  demonstrated  by  the  results  in  the  last  example,  the 
expected  convergence  rates  are  also  maintained  through  the  use  of  the  weight 


n 

m  =  0 

m  =  1 

m  =  2 

m  =  3 

It 

m  —  5 

m  =  6 

0 

61.96 

128.95 

211.56 

309.58 

422.56 

550.38 

692.75 

1 

241.23 

368.90 

513.02 

673.33 

849.82 

1042.07 

1249.92 

2 

540.46 

728.34 

932.93 

1154.26 

1392.14 

1646.34 

1916.70 

3 

959.46 

1207.39 

1472.15 

1753.95 

2052.63 

2367.90 

4 

1498.20 

1806.12 

2131.29 

2473.02 

5 

2156.69 

2524.45 

2909.31 

3311.57 

6 

2934.91 

3362.52 

3807.60 

4269.79 

Table  6.  Natural  frequencies  deriving  from  the  Bessel  expansions  (in  hertz). 


n 

m  =  0 

m  =  1 

m  =  2 

m  =  3 

3 

It 

m  =  5 

m  =  6 

0 

61.96 

128.95 

211.55 

309.52 

422.56 

550.38 

692.75 

1 

241.23 

368.96 

513.04 

673.40 

849.83 

1042.08 

1249.93 

2 

540.46 

728.35 

932.98 

1154.31 

1392.18 

1646.40 

1916.79 

3 

959.50 

1207.41 

1472.31 

1754.16 

2052.83 

2368.17 

4 

1498.37 

1806.28 

2131.37 

2473.57 

5 

2157.22 

2525.17 

2910.49 

3313.05 

6 

2936.35 

3364.46 

3810.17 

4273.26 

Table  7.  Natural  frequencies  obtained  via  the  Fourier-Galerkin  scheme  with  M  =  6,  N  =  24. 
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(m  =  0,  n  =  1) 


(m=  1,  n=  1) 


I 


(tn  =  2,  n  =  2)  (m  =  5,  n  =  0) 


Figure  4.  The  (0, 1),(1, 1),(2, 2)  and  (5,0)  modes  obtained  with  Bessel  expansions. 


(ni  =  0,  n=  1) 


(m- 


(tn  =  2,  n  =  2) 


(m  =  5,  D  =  0) 


Figure  5.  The  (0, 1),(2,2)  and  (5,0)  modes  obtained  via  the  Fourier-Galerkin  scheme 
with  Af  =  5  and  N  =  16. 
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4  Application  -  Structural  Acoustic  Noise  Control 

In  the  previous  two  sections,  Fourier-dalerkin  inethods  for  isolated  acoustic  and  stun  tural 
systems  have  been  developed  and  illustrated  through  a  s'-t  of  examples.  In  this  section,  a 
structural  acoustic  application  will  he  used  to  tlemonstrate  tlu*  maimer  though  whiih  these 
technicpies  can  he  comhined  to  yield  an  effective  approximation  scheme  for  simulating  the 
dynamics  of  more  complex  systems.  When  motleling  the  experimental  setup  which  motivates 
this  application,  care  was  taken  to  include  the  full  cou])ling  between  the  structural  dynamics 
and  the  enclosed  acoustic  field,  internal  damping  in  the  structure,  and  the  effects  of  actuators 
such  as  piezoceramic  patches  which  are  homled  to  the  structure.  The  full  inclusion  of  these 
components  leads  to  a  model  which  accurately  captures  the  physics  of  the  system  hut  is  diffi<  ult 
to  approximate  using  standard  modal  techni(iues  since  modes  for  the  coupled  system  (including 
coupling,  damping  and  actuator  effects)  must  themselves  be  approximated  since  they  are  not 
known  analytically.  I'sing  a  Fourier-Clalerkin  technique,  however,  the  approximate  system  can 
he  directly  constructed  with  several  of  the  components  actually  consisting  of  matrices  from  the 
isolated  acoustic  and  structural  systems.  lu  addition  to  illustrating  the  method  for  the  specific 
setup  described  here,  this  demonstrates  the  flexibility  of  the  method  for  approximating  the 
dynamics  for  general  structural  acoustic  systems  having  circular  or  cylindrical  geometries. 


4.1  The  Structural  Acoustic  System 

The  structural  acoustic  system  de.scribed  here  models  an  experimental  setup  currently  being 
used  for  validation  experiment.s  in  the  Arunstirs  Division,  NASA  Langley  Research  Outer,  and 
the  geometry  and  physical  specifications  were  chosen  so  as  to  be  consistent  with  that  apparatus. 
Specifically,  the  experimental  apparatus  is  modeled  by  a  cylindrical  domain  fl  having  length  f 
and  radius  a  as  pictured  in  Figure  6.  At  one  end  of  the  cylinder  is  a  clamped  flexible  plate  of 
thickness  k  which  is  assumed  to  have  Kelvin- Voigt  damping.  Bonded  to  the  plate  are  sectorial 
piezoceramic  patches  which  are  placed  in  pairs  and  excited  out-of-phase  so  that  a  bending 
moment  is  produced  when  voltage  is  appled.  The  patches  and  glue  layer  are  cussumed  to  have 
thicknesses  T  and  Tw,  respectively. 


Figure  6.  (a)  The  cylindrical  acoustic  cavity;  (b)  The  circular  plate  with  patches. 
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As  discussed  in  [3,  4],  an  appropriate  linearized  variational  ((Mier^y)  form  id  tlic  (Dupled 
system  equations  modeling  this  setup  is 


/  +  /  p/V(j!)  ■  4-  f  g^du; 

Jo  Jo  Jo 


{4.1«) 


+  l^p,Uw.,T,dT  +  (I.H) 

+  /  Pf  {4>tV  -  wtt)  d'f  (  l.lr) 

^Vq 

=  f  Y^l^fui{t)Xx{r,e)\/h]d'y-\-l  fgd')  (4. Id) 

•'To  •'To 

for  all  appropriate  test  functions  ^  and  t]  (this  will  be  clarified  below).  The  differentials  here 
are  du)  =  rdrdOdz  and  d7  =  rdrdO  with  the  overbars  again  denoting  complex  conjugation. 
We  now  consider  the  components  of  the  equations  and  compare  with  those  found  previously. 
Throughout  this  section,  the  subscripts  p  and  c  will  be  used  to  denote  plate  and  wave  (cavity) 
components,  respectively. 

Acoustic:  With  €  L^{Q)  again  denoting  the  acoustic  velocity  potential,  the  wave  dynamics 
are  contained  in  (4. la)  where  the  test  functions  ^  are  elements  in  H\fl)  (see  Section  2 
for  details  concerning  the  Hilbert  spaces  for  this  problem).  When  approximating  the 
dynamics  of  the  system,  the  2-D  basis  described  in  Section  2  is  tensored  with  an  axial 
Legendre  basis  to  yield  the  approximate  solution 

Pc  Me 

4.^{t,r,e,z)  =  Yl  E  E  . 

p=:0  7U=  — Mr  n=0 

P+|m|+n^O 


where  again. 


rn 


r  m 

~l  •'i 


jm!  =  o,---,r> 

|m|  =  6,  •  •  •,  A/c 


and 


f  A(»‘)-l/3  7  p  =  m  =  0,n=l  f  jVr+l  ,  p+(m|^0 

^^’”‘(0  =  i  7  Nr^  =  { 

I  ^n(^)  7  otherwise  (  ,  p  =  7/1  =  0 

Here  Pn(j’)  and  Pp{z)  are  the  and  p‘^*  Legendre  polynomials  which  have  been  mapped  to 
the  intervals  (0,  a)  and  (0,^),  respectively.  We  remind  the  reader  that  the  term  Pi{r)—  1/3 
when  p  =  m  =  0,n  =  1  results  from  enforcing  the  condition  /q  (i)-^{t.,r,9,  z)du}  =  0  so  as 
to  guarantee  that  the  functions  are  suitable  as  a  basis  for  the  quotient  space  while  /  '’'‘I 
incorporates  the  analytic  decay  at  the  origin  with  m  truncated  to  control  conditioning. 
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Structural:  Equation  (4.1b)  cuntaiiis  the  internal  plate  dynainies  while  the  external  cuiitribu- 
tions  due  to  the  pertvirbing  noise  source  /  and  the  t*xritation  of  the  j)atclies  are  contained 
in  (4.1(1)  (compare  with  (-4.11)  in  the  last  section).  As  was  the  case  when  th<‘  isolatcsl 
plate  was  considered,  iv  G  L^lVo)  denot(*.s  tlie  transvt^rse  displacement  and  the  test  tunc- 
tions  ;/  are  taken  in  A/o(ro).  When  discretizing  the  system,  the  plate  displacement  is 
approximated  by  tL'-^{t.,r,6)  as  defined  in  (4.12)  with  A/^  and  used  lu're  to  denote  the 
Fourier  and  spline  limits,  respectively. 

Coupling:  The  coupling  between  the  plate  dynamics  and  interior  acoustic  field  is  incorporated 
by  including  the  backpressur*'  pj<j)i(t^i\0)  from  the  interior  field  as  a  force  on  the  plate 
and  assuming  the  continuity  of  velocity  condition 

=  —  Wt{t,r.O) 

()z 

at  the  surface  of  the  plate  (this  latter  condition  is  often  designated  the  momentum  con¬ 
dition  in  the  literature).  In  the  weak  form  of  the  .system  equations,  these  conditions  are 
manifested  as  the  first  and  second  terms,  respectively,  in  (4.1c). 

With  the  components  thus  described,  we  are  now  in  a  position  to  form  the  matrix  system 
which  results  when  the  system  dynamics  are  approximated.  The  Af  =  (2 Afp  +  1 )( A^p  +  1 )  —  1  and 
M  =  {2Mc  -f  1)(A^(.  +  l)(/^f  +  1)  —  I  dimensional  approximating  plate  and  cavity  subspaces  are 
taken  to  be  =  spanlB,^}^,  and  =  spanlBA*}^,  ,  respectively,  where  Bf  and  are 
the  plate  and  cavity  bases  described  above.  Defining  V  =  Af  +  AA,  the  approximating  state 
space  is  x  and  the  product  space  for  the  first  order  system  is  JC  —  x  . 

The  restriction  of  the  infinite-dimensional  system  (4.1)  to  x  then  yields  the  matrix 
system 

=  A'^y'^(t)  +  B'^u{t)  +  F’’(<) 

(4.2) 

M'^y'^(O)  =  y^  . 

Here  y^{t)  —  [d^{t),  W^(t),  d^(t),  VV^(t)]^,  with  the  components  i?^(<)  =  [(j>^ {t) ^ 

and  W^(t)  =  [u)f(t),  •  •  • ,  u;]5)-(t)]^,  denotes  the  approximate  state  vector  coefficients  while 

u{t)  =  [ui(t),  •  ■  • ,  Us(t)]^  contains  the  s  patch  input  variables.  The  system  matrices  and  vectors 


M^  = 


A^  = 


^t9z 


-n.p  ^c2 


=  \  0  0 


OBf 


F'^(t)  =  [  0  0  Gf  (0  Ffit) 


We  now  separate  the  matrices  into  tho.se  containing  acoustic,  structural  and  coupling  compo¬ 
nents  to  better  illustrate  the  connections  with  those  developed  for  the  isolated  acoustic  and 
structural  problems  in  the  previous  sections. 
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Structural:  The  Af  x  ,V  component  matrices  A'^j  aiul  A';'^  are  the  niass.  stiffness  and 
damping  matrices  which  arise  when  solving  the  dampeil  plate  ecpiation  with  fixed  hound* 
ary  conditions  while  the  x  1  vectors  and  Ff{f)  are  the  corresponditig  control  and 
forcing  terms  (see  (3.13)  for  the  various  definitions). 

Acoustic:  The  AA  x  AA  matrices  and  are  the  mass  and  stiffness  matrices  which 

arise  when  solving  the  uncoupled  wave  equation  with  Neumann  houndarv  conditions  on 
a  cylindrical  domain.  As  a  results  of  the  tensor  propertie*s  of  the  3-1)  l)asis.  they  ran  lie 
succinctly  defined  as  the  tensor  products 

w  , 

where  and  are  given  in  (2.15)  and  the  (P^-  -f  1)  x  (P,.  4-  1)  matrices  .\P  and  A'_. 
are  defined  by 

=  /  Pd-)PA-)d^  ,  =  j('  a'(--)A;(--R- . 

The  construction  of  and  is  completed  by  updating  the  row  and  column  affected 
by  the  alterations  used  to  guarantee  that  the  functions  are  a  basis  for  the  quotient  space. 
Finally,  G^{t)  is  a  forcing  term  which  incorporates  any  acoustic  sources  (see  (2.15)). 

Coupling:  The  contributions  from  the  coupling  terms  are  contained  in  the  matrices 

where  Bf  and  B^  are  plate  and  acoustic  basis  functions,  respectively,  and  the  index 
ranges  are  k,i  =  1,  •  •  • ,  A1  and  r,p  =  T  •  ’  •  j-Af. 

Initial  Conditions:  The  vector  , gf , (see  (2.15)  and  (3.13))  contains  the 

projections  of  the  initial  values  into  the  approximating  finite  dimensional  subspaces. 

We  point  out  that  several  of  the  component  matrices  are  identical  to  those  defined  when 
considering  the  isolated  acoustic  and  structural  systems  while  others  can  be  efficiently  con¬ 
structed  from  those  components  through  basic  operations  such  as  tensor  products.  In  fact, 
only  the  coupling  matrices  and  containing  quadrature  values  for  the  plate  and  acoustic 
bases  must  be  constructed  solely  for  this  problem.  This  potential  for  decomposition  into  exist¬ 
ing  structural,  acoustic  and  coupling  matrices  also  exists  for  more  complex  structural  acoustic 
systems  when  a  Fourier-Galerkin  approximation  scheme  is  used  and  is  a  further  advantage  of 
the  method  in  coupled  problems  involving  linear  or  weakly  nonlinear  components. 

4.2  Example  5 

As  a  first  test  demonstrating  the  convergence  of  the  method  for  the  fully  coupled  system,  we 
considered  the  problem  (4.1)  with  the  dimensions  a  =  .6,  h  =  .00127,  =  1.1  and  parameter 

choices  pf  =  l.21,c  —  343, pp  =  2700,*/  =  .33,0  =  =  13.6007  and  cp  =  .00011222  (the 
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effects  of  the  patches  are  iieglec  led  in  this  example  hut  will  be  incorporated  in  Fxample  (i). 
The  true  solutions 

w{r,9)  =  t:\coH('>wr/a)  —  l)sin^ 

/■,  0,  z)  =  —  ^/^(cos(27rr/«)  —  \  )z[(  —  z)^  sin  0 

were  used  to  generate  forcing  functions  /  and  </  on  the  plate  and  in  the  cavity,  respe(  lively  (see 
the  convergence  tests  in  .Sections  2  and  The  resulting  matrix  system  (4.2)  was  then  integrate*d 
to  time  T  =  A  and  the  absolute  and  relative  errors  in  the  plate  displacement,  potential  ainl 
pressure  at  that  time  were  calculated.  The  results  for  four  sets  of  Clalerkin  limits  are  reported 
in  Tables  8,  9  and  10.  In  these  tables,  the  subscripts  p  and  c  in  the  tablets  again  refer  to  the 
plate  and  wave  (cavity)  indices,  respectively,  and  the  measurements  were  made  on  a  80  x  80 
grid  on  the  plate  and  a  10  x  10  x  10  grid  in  the  cavity.  As  demonstrated  by  the.se  results,  the 
method  for  the  fully  coupled  system  is  exhibiting  a  convergence  rate  simi'ar  to  that  noteii  when 
the  individual  components  were  tested.  Moreover,  sufficient  accuracy  is  obtained  with  tractable 
matrix  sizes  so  that  physical  details  are  captured  when  the  approximation  method  is  used  to 
simulate  system  dynamics  and  develop  feedback  control  techniques.  Finally,  we  point  out  that 
in  spite  of  the  fact  that  we  are  approximating  nonaxisymmetric  solutions,  the  computations  for 
this  example  could  be  performed  on  a  -Sparc  1 0-class  workstation. 


Mp  Np  Me  Ne  Pc  size(Af^) 

1  5  i  2  2  86  X  86  .4366  -  5  .2207  -  2 

1  5  1  4  4  182  x  182  .4294-  5  .2171  -2 

1  10  1  6  6  356  x  356  .2674-  6  .1352-3 

1  20  1  8  8  608x  608  .1588-  7  .8030-5 


Table  8.  Absolute  and  Relative  Errors  in  Plate  Displacement  at  T  =  .1. 


A/p  Ap  Me  Ne  Pc  sizeiA'^)  \\<t>true  -  <t>apv\\ 

"1  5  i  2  2  86  X  86  .3132-  2  .3569  -  0 

1  5  14  4  182  x  182  .1545-3  .1760-1 

1  10  1  6  6  356x  356  .9554-5.  .1089-2 

1  20  1  8  8  608  x  608  .2251  -6  .2.565-4 


Table  9.  Absolute  and  Relative  Errors  in  Potential  at  T  =  .1. 


A/p  Ap  Me  Ne  Pe  sizc(A^)  \\ptrue  -  PappW 

"1  5  i  2  2  86x  86  .7573  -  1  ^3.566  -  0 

1  5  1  4  4  182x  182  .3739-2  .1761  -1 

1  10  1  6  6  3.56x  3.56  .2312-3  .1089-2 

1  20  1  8  8  608  X  608  .5453  -  5  .2568  -  4 


Table  10.  Absolute  and  Relative  Errors  in  Pressure  at  T  =  .1. 
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4.3  Example  6 

In  this  final  example,  the  numerical  techniciues  just  descrihcd  are  used  to  simulate  the  system 
dynamics  of  the  experimental  structure  described  at  the  beginning  of  this  section.  A  common 
technique  for  determining  the  dynamics  and  natural  fretiuencie,-.  of  coupled  structural  acoustic 
systems  is  to  subject  the  structure  to  an  impulse  (through  an  impact  hammer  hit  or  a  voltage 
spike  to  the  patch)  and  measure  the  resulting  tim*'  and  freriuency  r(‘spon.ses.  By  ('xciting  the 
system  in  this  manner,  a  wide  spectrum  ot  frec|uen<ies  cati  be  excited  and  vaiious  systtun 
properties  determined.  Because  there  is  no  acoustic  souix'e  inside  the  cavity.  </  =  0  iti  {  1.1a). 

The  hammer  impact  can  be  numerically  simulated  by  a  triangular  force  f  to  a  |)oint  on  tin¬ 
plate  while  the  voltage  spike  can  be  approximated  by  a  short  iluration  triangular  vcjltage  u:  in 
both  ca.ses,  tiie  approximate  system  response  is  calculated  via  (4.2). 

In  order  to  remain  consistent  with  the  experimental  setup  Ix-ing  modeled,  tin-  h-ngtli  and 
radius  of  the  cavity  were  taken  to  be  1.0fi68m  (42")  and  a  =  .2286 //(  (9").  n-spectivel>  with  a 
plate  having  thickness  h  =  .00127  m  (.05")  mounted  at  one  end.  A  pair  of  circular  piczoceramic 
patches  having  thickne.ss  T  =  .0001778m  (.007")  and  radius  rad  =  .OlOOorn  (.75")  were  located 
at  the  center  of  the  plate  (see  fdgure  7).  The  physical  parameters  that  were  chosen  for  tin- 
structure  and  acoustic  cavity  are  summarized  in  Table  1 1.  The  flexural  rigidity  D  for  the  plate 
was  obtained  using  the  “handbook”  value  E  =  7. 1  x  lO’^.'V/m.^  for  the  Young’s  modulus  of 
aluminum.  The  remaining  choices  are  comparable  to  values  found  when  estimating  paramete-rs 
for  the  isolated  plate  with  a  similar  patch  configuration  [2].  We  re-emphasize  that  in  general, 
these  parameters  must  be  determined  through  parameter  estimation  technic|ues  in  order  for 
the  PDE  model  to  fit  the  actual  physical  system  under  consideration.  As  demonstrated  in  [2], 
the  Fourier-Clalerkin  scheme  presented  here  performed  well  when  incorporated  in  a  fit-to-data 
routine  for  determining  physical  plate  parameters  thus  yielding  a  model  which  could  be  used 
for  implementing  model-based  feedback  control  techniques. 


Figure  7.  The  acoustic  cavity  with  a  pair  of  centered  circular  patches,  the  impact  point 
71  =  (.13,  ^)  and  the  observation  points  pi  =  (.125,  ^)  and  cl  =  (.125,  .05). 
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Structure 


Struc 

ture 

Acoustic  Cavity 

Parameter 

Plate  Plate  -H  Pzt 

Parameter 

Cavity 

Pp  ■  Tliicknesti  {kg/iiP) 

3.429  3.4S9 

Pj  (A'i//m'*) 

1.21 

D  {N  -m) 

C[)  (  N  •  7)1  •  .see) 

n 

fC^  {N/V) 

13.601  13.901 

1.150-4  2.250-4 

.33  .32 

.0267 

c  (m/sec) 

343 

Table  11.  Physical  Parameters  for  the  structure  and  acoustic  cavity. 


Voltage  Spike  to  a  Centered  Patch  Pair  -  Axisymmetric  Response 

As  a  first  example  demonstrating  the  response  of  the  system  model  to  a  simulated  impact, 
we  applied  a  triangular  voltage  spike  u{t)  as  input  to  the  patches  and  integrated  the  matrix 
system  to  obtain  a  time  history  of  the  response.  Because  the  patches  are  centered  on  the 
plate,  the  response  was  uniform  in  6  which  implied  that  the  Fourier  limits  could  be  taken  to 
be  Mp  —  Me  =  0.  The  remaining  limits  were  taken  to  be  Np  =  12,  =  9  and  Pc  =  9 

(these  values  resolve  the  range  of  frequencies  being  examined).  A  time  history  of  the  system 
response  throughout  the  interval  [0,  .5]  was  calculated  at  the  plate  and  cavity  points  pi  = 
(.125,  and  cl  =  (.125,  ^,.05)  (see  Figure  7),  and  the  resulting  trajectories  and  frequencies 
are  plotted  in  Figure  8.  Tliis  temporal  interval  was  chosen  since  it  was  sufficiently  long  so  as 
to  demonstrate  the  system  dynamics  but  short  enough  so  that  the  higher  frequency  responses 
were  not  completely  lost.  The  off-center  observation  points  were  chosen  to  demonstrate  the 
generality  of  the  method  and  provide  a  basis  for  comparison  with  the  nonaxisymmetric  resc.its 
in  the  next  example.  Finally,  the  system  frequencies  are  summarized  in  Table  12  with  the 
notation  p  and  c  being  used  to  designate  those  frequencies  which  are  observed  at  the  plate 
point  pi  and  cavity  point  cl,  respectively. 

To  determine  the  effects  of  coupling,  internal  plate  damping,  and  the  presence  of  the  patches 
on  the  system,  it  is  illustrative  to  compare  the  system  results  listed  in  Table  12  with  those  of 
the  isolated  components.  In  Section  3,  the  natural  frequencies  for  an  isolated  and  undamped 
plate  were  obtained.  The  natural  frequencies  for  the  isolated  wave  equation  in  a  cylindrical 
cavity  can  be  determined  by  separating  variables  in  the  3-D  wave  equation  having  Neumann 
boundary  conditions.  This  leads  to  a  3-D  Helmholz  equation  which,  after  an  analysis  similar 
to  that  presented  in  Section  2,  yields  the  natural  frequencies 


where  p  =  1,2,  •••  ,  t?i  =  0,  1, 2,  •  •  •  ,  n  =  0,1,2,  •••  and  A„„i  =  are  zeros  of  (2.6)  (see 

[3,  6,  9]  for  details).  For  the  previously  mentioned  dimensions,  the  frequencies  of  axisymmetric 
modes  (m  =  0)  are  listed  along  with  those  for  the  isolated  and  undamped  plate  in  Table  13 
(see  Tables  2  and  6  for  a  compilation  of  the  zeros  7„i,i  and  natural  frequencies  for  the  isolated 
plate,  respectively). 


By  comparing  the  values  for  the  individual  plate  and  cavity  in  Table  Id,  with  the  system 
harmonics  in  Table  12,  it  can  be  seen  that  although  the  frt'ciuencies  agree  (|uite  closely,  then* 
are  slight  differences  due  to  the  fact  that  the  system  (4.1)  involves  not  only  cou|ding  between 
the  plate  and  cavity  but  also  includes  damping  in  the  jilate.  .Specitically,  the  system  fie(|uen- 
cies  associated  with  strong  plate  responses  tend  to  be  slightly  less  than  those  of  the  uncoupled 
and  undamped  plate  while  those  system  fre(|uencies  as.sociatetl  with  strong  acoustic  responses 
are  slightly  higher  than  the  natural  acoustic  frequencies  of  an  isolated  hardvvalled  cavity  hav¬ 
ing  these  dimensions.  Hence,  through  the  Fourier-CJalerkin  method,  physical  characteristics 
about  the  system  are  determined  and  differences  between  the  couphHl  system  res|)onse  which 
includes  damping  and  the  isolated,  undamped  components  are  illustrated.  This  latter  obser¬ 
vation  also  indicates  the  difficulties  which  would  be  encountered  if  one  were  to  attem|jt  to 
simulated  the  coupled  system  dynamics  using  modes  for  the  isolated  systems. 


Natural  .System  Frequencies 


59 

P>f 

164 

c 

915 

p 

240 

p.c 

324 

c 

929 

p,<' 

537 

P.f 

483 

c 

971 

c 

645 

c 

807 

Table  12.  System  frequencies  obtained  with  Mp  =  0,  A^p  =  \2,Mc  =  0,  A^o  =  and  Pp  =  9 
basis  functions;  p  -  frequencies  observed  at  the  plate  point  pi  =  (.125,47r/;i), 
c  -  frequencies  observed  at  the  cavity  point  cl  =  (.125,47r/3,  .05).  System  freq¬ 
uencies  can  be  compared  with  the  ordered  frequencies  of  the  isolated  and  undamped 
plate  and  cavity  given  in  Table  13. 


Plate  (/,„„) 

Wave 

.  fmnp  ) 

(0,0) 

62 

(0,0,1) 

161 

(0,1,0) 

915 

(0,1) 

241 

(0,0,2) 

322 

(0,1,1) 

929 

(0,2) 

540 

(0,0,3) 

482 

(0,1,2) 

9.59 

(0,0,4) 

643 

(0,0,5) 

804 

(0,0,6) 

965 

Table  13.  Axisymmetric  natural  frequencies  for  the  isolated  and  undamped  plate  and  cavity 
(in  hertz). 
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Figure  8.  The  plate  and  pressure  responses  to  a  voltage  spike. 


Noncentered  Hammer  Impact  -  Non-axisymmetric  Response 

As  indicated  previously,  a  second  means  of  exciting  the  system  is  through  a  hammer  impact 
which  can  be  modeled  by  a  short  triangular  input  force.  To  demonstrate  the  approxima¬ 
tion  of  general  nonaxisymmetric  system  dynamics,  this  force  impulse  was  applied  at  the  point 
71  =  (.130,  |),  and  a  time  history  of  the  system  response  for  the  time  interval  [0,  .5]  was  cal¬ 
culated  at  the  plate  and  cavity  points  pi  =  (.125,  and  cl  =  (.125,  .05),  respectively  (see 

Figure  7).  The  resulting  acoustic  pressure  and  plate  acceleration  are  plotted  along  with  the  cor¬ 
responding  frequency  responses  in  Figure  9  (the  plate  acceleration  models  data  that  would  be 
experimentally  obtained  with  an  accelerometer,  [2],  and  more  clearly  demonstrates  the  higher 
frequencies  than  does  the  plate  displacement).  To  resolve  the  frequencies  below  1000  hertz, 
the  basis  limits  were  taken  to  be  Mp  =  4,  Ap  =  12,  Me  =  2,  =  9  and  Pp  =  9  which  yielded 

1230  coefficients  in  the  system  (4.2).  As  demonstrated  by  this  choice  of  limits,  the  number  of 
Fourier  coefficients  in  the  plate  and  wave  expansions  can  be  chosen  to  differ  in  order  to  reduce 
system  sizes  although  care  must  be  taken  when  doing  .so  since  the  system  response  contains 
contributions  from  both  the  plate  and  wave  components. 

The  system  frequencies  are  summarized  in  Table  14  where  again,  the  notation  p  and  c  are 
used  to  designate  those  frequencies  which  are  observed  at  the  plate  and  cavity  point  pi  and  cl, 
respectively.  For  comparison,  corresponding  natural  frequencies  for  an  isolated  and  undamped 
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plate  (see  Table  6)  and  an  isolated  acoustir  cavity  (see  (4..'}))  are  compiled  in  ral)le  lo.  By 
comparing  the  results  in  these  table,  it  can  be  notetl  that  six  system  modes  ( orrespondiiig  to 
those  of  the  isolated  plate  are  excited  with  four  of  these  responses  being  observed  at  l)otli  the 
plate  and  cavity  points  (we  emphasize  that  because  these  are  truly  system  modes,  the  remaining 
frequencies  of  240  hertz  and  422  hertz  can  also  be  measured  in  the  cavity  at  various  points:  tlie 
response  is  simply  weak  at  the  point  cl).  Similarly,  nineteen  system  modes  corresponding  to 
those  of  the  isolated  acoustic  cavity  are  excited  by  the  nonsymmetric  impact  with  the  strongest 
response  in  the  asymmetric  system  mode  having  a  frequency  of  460  hertz  (this  corresponds  to 
the  (1,0,0)  mode  for  the  isolated  cavity).  As  was  noted  in  the  previous  example  demonstrating 
the  symmetric  excitation,  the  system  frequencies  associated  with  strong  plat<'  responses  t<'nd 
to  be  slightly  less  than  those  of  the  uncoupled  and  undamped  plate  while  those  system  fre¬ 
quencies  associated  with  strong  acoustic  responses  are  slightly  higher  than  the  natural  acoustic 
frequencies  of  an  isolated  hardwalled  cavity  having  commensurate  dimensions.  We  reiterate 
that  these  differences  are  due  to  the  coupling  between  the  plate  and  cavity  as  well  as  damping 
in  the  plate.  This  again  provides  motivation  for  considering  a  general  Fourier-Clalerkin  method 
of  the  type  described  here  when  approximating  the  dynamics  of  a  complex  coupled  system. 
The  effect  of  coupling  and  damping  are  implicitly  included  in  this  approach  while  they  must  be 
explicitly  (and  occasionally  artificially)  enforced  when  employing  modal  approaches  using  the 
eigenfunctions  of  the  isolated  components. 


Natural  Frequencies  of  the  Coupled  Structural  Acoustic  System 


Table  14.  System  frequencies  obtained  with  Mp  =  4,  A^p  =  12,  =  2,  Nc  =  9  and  Pc  =  9 

basis  functions;  p  -  frequencies  observed  at  the  plate  point  pi  =  (.125,47r/d), 
c  -  frequencies  observed  at  the  cavity  point  cl  =  (.125,4^/3,  .05).  System  fre¬ 
quencies  can  be  compared  with  the  ordered  frequencies  of  the  isolated  and  un¬ 
damped  plate  and  cavity  given  in  Table  15. 


Plate 

{fvitl  ) 

Wavp  (/iuTip) 

62 

212 

(0,0,1) 

161 

||ran| 

929 

(1,0,5) 

916 

■ran 

241 

(2,1) 

513 

(0,0,2) 

322 

1  Q 

970 

(2,0,0) 

729 

540 

(2,1) 

933 

(0,0,3) 

482 

1 

439 

(2,0,1) 

750 

■mBI 

959 

310 

I  1 

643 

1  |ra 

468 

797 

(1,0) 

129 

(3,1) 

673 

1  qQb  I 

804 

1 

545 

874 

(1,1) 

369 

423 

1  ^0  1 

965 

1  Q 

653 

Iwpiill 

972 

(1,2) 

728 

(4,1) 

850 

1  nn  1 

915 

IliiaE 

779 

Table  15.  Nonaxisymmetric  natural  frequencies  for  the  isolated  and  undamped  plate  and 
cavity  (in  hertz). 
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Pkat«  Acc«tora«on  at  |  125  4pi/3t 


Amplitucte  Spactrum  ol  Dtsplac«mani  at  i  i2i>  4pi/3f 


Figure  9.  The  plate  and  pressure  responses  to  an  impart  at  (.13, 7r/2). 

5  Concluding  Remarks  and  Control  Applications 

In  this  paper,  a  Clalerkin  method  for  api^roximating  the  dynamics  of  linear  or  weakly  nonlinear 
systems  in  circular  geometries  is  presented.  While  the  method  is  presented  for  linear  acoustic, 
structural  and  structural  acoustic  problems,  the  technic|ue,s  are  sufficiently'  general  so  as  to  be 
applicable  in  a  large  variety  of  applications  involving  the  approximation  of  linear  or  weakly 
nonlinear  operators  on  circular  domains. 

To  take  advantage  of  periodicity,  Fourier  expansions  are  first  used  to  approximate  the  tan¬ 
gential  behavior  of  the  solution.  The  crux  of  the  method  then  involves  the  construction  of 
radial  basis  functions  by  combining  traditional  approximating  functions  (finite  elements,  piece- 
wise  splines,  spectral  polynomials)  with  terms  of  the  form  where  ih  is  a  truncated  Fourier 
coefficient  and  the  coordinate  singularity  is  assumed  to  be  at  the  origin.  This  latter  term  incor¬ 
porates  the  decay  properties  of  the  Bessel  or  analytic  solution  near  the  origin  and  is  included 
so  as  to  guarantee  that  the  approximate  .solution  is  unique  and  differentiable  at  the  coordi¬ 
nate  singularity.  The  truncation  of  the  Fourier  coeHicient  pnwides  a  means  of  ct)ntrolling  the 
conditioning  of  the  resulting  matrix  system. 

As  demonstrated  by  several  numerical  examples,  the  resulting  Fourier-dalerkin  technique 
is  accurate,  efficient  and  sufficiently  flexible  so  as  to  be  easily  extended  to  complex  coupled 
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systems.  By  varying  the  rhoire  of  radial  basis  fuiKtioiis.  the  method  can  lx-  tailored  to  a 
variety  of  applications  and  boundary  conditions  while  the  inclusion  of  the  dec  ay  c  cuiipoueiit  /  '  • 
guarantees  that  accuracy  is  maintained  throughout  the  entire  circular  domain  for  sullic  icmt  ly 
smooth  system  inputs.  Hence  the  Legendre’  poivnomials  atnl  moditied  cidiic  splines,  iisc’d  in 
this  presentation  to  approximate  acoustic  and  structural  dynamics,  can  easily  be  replac  es!  Icy 
other  spectral  functions,  splines  or  finite  ehunents  as  warranted  by  the-  sit'iation. 

A  final  advantage  of  the  method  arises  when  the  Fourier-dah'rkin  teclmicjiies  are  uscsl  to 
compute  feedback  gains  in  optimal  control  problems.  As  discus.sed  in  [5],  where’  the-  me'tluKl  was 
incorporated  in  a  control  scheme  for  re*ducing  noi.se  in  a  cylindrical  couplexl  struc  tural  acoustic- 
system,  the  optimal  control  problem  for  reducing  a  periodic  ally  driven  state  c  in  the  state*  sjcac  e’ 
'H  can  be  posed  under  approximation  as  that  of  determining  a  suitable  u  in  a  control  sjeae  e-  I 
which  minimizes 

JV)  =  ^  r^.l) 

subject  to  satisfying  the  matrix  system  of  dimension  V  which  results  when  the  e(|uations 
modeling  the  system  dynamics  are  discretized.  Here  r  is  the  period  anel  and  are  matrices 
which  can  be  used  to  weight  various  components  of  the  approximate  state  and  control.  As  noted 
in  [5],  a  suitable  choice  for  ,  when  an  energy  inner  product  is  associated  with  the  state  spac  e 
'H,  is  a  diagonal  multiple  of  the  mass  matrix  which  results  when  a  (lalerkin  method  is  used 
to  discretize  the  weak  form  of  the  system  equations.  Since  the  mass  matrix  c  onstruc  ted  in  this 
manner  is  the  identity  with  respect  to  the  energy  inner  product,  the  choice  =  VM^ ,  V  a 

diagonal  weighting  matrix,  in  (5.1)  minimizes  a  weighted  measure  of  the  state  energy.  Hence  a 
Fourier-Galerkin  method  of  the  type  described  here  can  also  be  advantageous  when  calculating 
feedback  gains  for  experimental  and  numerical  implementation  of  optimal  control  techniques 
to  problems  posed  on  circular  geometries. 
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